Determination of the stiffness of the nuclear symmetry energy from isospin diffusion 
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With an isospin- and momentum-dependent transport model, we find that the degree of isospin 
diffusion in heavy ion collisions at intermediate energies is affected by both the stiffness of the nuclear 
symmetry energy and the momentum dependence of the nucleon potential. Using a momentum 
dependence derived from the Gogny effective interaction, recent experimental data from NSCL/MSU 
on isospin diffusion are shown to be consistent with a nuclear symmetry energy given by i5sym(p) ~ 
31.6{p/ poY '^^ at subnormal densities. This leads to a significantly constrained value of about —550 
MeV for the isospin-dependent part of the isobaric incompressibility of isospin asymmetric nuclear 
matter. 
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Knowledge on the density dependence of nuclear sym- 
metry energy is important for understanding not only the 
structure of radioactive nuclei P, S IS 0L but also many 
important issues in nuclear astrophysics [^Q, such as nu- 
cleosynthesis during pre-supernova evolution of massive 
stars and the cooling of protoneutron stars. Although 
the nuclear symmetry energy at normal nuclear matter 
density is known to be around 30 MeV from the empirical 
liquid-drop mass formul a M , M , its values at other den- 
sities are poorly known j^, lnj ■ Studies based on either 
microscopic many-body theories or phenomenological ap- 
proaches have so far given widely divergent predictions on 
the density dependence of nuclear symmetry energy . 
Empirically, the incompressibility of asymmetric nuclear 
matter is essentially undetermined [l2l |. even though the 
incompressibility of symmetric nuclear matter at its satu- 
ration density po = 0.16 fm^'' has been determined to be 
231 ±5 MeV from nuclear giant monopole resonances [l^ 
and the equation of state at densities of 2po < p < 5po 
has been constrained by measurements of collective flows 
in nucleus- nucleus collisions T^. 

In light of the new opportunities provided by radioac- 
tive beams, a lot of interests and activities have recently 
been devoted to extract information on the density de- 
pendence of nuclear symmetry energy from reactions in- 
duced by such nuclei [HEElElElliSIIllillil 
I24, 25, 26]. In particular, isospin diffusion in heavy- 
ion collisions is found to depend sensitively on the den- 
sity dependence of nuclear symmetry energy [27I |2^ . 
Within a momentum-independent transport model, in 
which the nuclear potential depends only on local nu- 
clear density, the isospin diffusion data from recent ex- 
periments at the NSCL/MSU (National Superconducting 
Cyclotron Laboratory at Michigan State University) was 
found to favor a quadratic density dependence for the 
interaction part of nuclear symmetry energy iflj. This 
conclusion has stimulated much interest because of its 
implications to nuclear many-body theories and nuclear 



astrophysics. However, the nuclear potential acting on a 
nucleon is known to depend also on its momentum. For 
nuclear isoscalar potential, its momentum dependence is 
well-known and is important in extracting information 
about the equation of state of symmetric nuclear matter 
[IlllllllllSlllllllllllllli. Very recently, the 
momentum -dep endence of the isovector (symmetry) po- 
tential 0, 139I I40I liH was also shown to be important 
for understanding a number of i sosp in related phenomena 
in heavy-ion reactions QS l44j . It is thus necessary 
to include momentum dependence in both the isoscalar 
and isovector potentials for studying the effect of nuclear 
symmetry energy on isospin diffusion. In this Letter, we 
shall show that the isospin diffusion data are consistent 
instead with a softer symmetry energy that is nearly lin- 
ear in density within the momentum-dependent effective 
interaction used in present study. 

Our study is based on an isospin-dependent transport 
model IBUU04 that uses experimental nucleon-nucleon 
cross sections in free space and includes the momen- 
tum dependence in both the isoscalar and isovector po- 
tentials Although the momentum dependence of 
the isoscalar potential is known empirically, that of the 
isovector potential is not as well determined In the 
IBUU04 model, it is based on the Gogny effective interac- 
tions. Specifically, the potential U{p, 6, p, r) for a nucleon 
with isospin r (1/2 for neutrons and —1/2 for protons) 
and momentum p in asymmetric nuclear matter at total 
density p is given by 
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where pr and p'^ denote proton or neutron density with 
T ^ t'; and 5 = (p„ — Pp)lp is the isospin asymmetry. 
The fr[v,-p) denotes the phase-space distribution func- 
tion at coordinate r and momentum p. The correspond- 
ing momentum-dependent interaction (MDI) leads to an 
incompressibihty of A'o = 211 MeV for the symmetric 
nuclear matter at saturation density. On the right hand 
side of Eq.(P), the first four terms with a = 3/4 and 
B — 106.35 MeV describe the momentum-independent 
interaction. The terms with parameters Cr,T — —11.7 
MeV and Cr.r' = —103.4 MeV describe the momentum- 
dependent interaction of a nucleon of isospin r and mo- 
mentum p with like and unlike nucleons in the back- 
ground fields, respectively. With the parameter A = 
l.Op^p, where p^p denotes nucleon Fermi momentum at 
po, the isoscalar potential {Un{p,p) + Up{p,p))/2 coin- 
cides with predictions from the variational many-body 
theory using inputs constrained by nucleon-nucleon scat- 
tering data j4^, and the isovector potential (Un{p,p) — 
Up{p,p))/2 also agrees with the momentum dependence 
of the Lane potential extracted from low-energy nucleon- 
nucleus scattering experiments |42l | and recent neutron- 
nucleus data at 96 MeV H. 



80 



60 



> 

(U 

^40 



20 





/ 


MDI interaction 






/ 




/ 








^ • ■ - ' ' 












x^T" — ■ 










/y y 




/y ■ ' . ^ ' 









0.0 



0.5 



1.0 

p/Po 



1.5 



2.0 



FIG. 1: (Color online) Density dependence of nuclear sym- 
metry energy for different values of x parameter in Eq.Q. 



The parameter x in Eq.(Q is introduced to mimic the 
various theoretical predictions for the density dependence 
of the nuclear symmetry energy i^symCp), which is defined 
via the parabolic approximation to the nucleon specific 
energy in an isospin asymmetric nuclear matter [3, 0| , 
i.e., 



E{p, 8) = E{p, ,5 - 0) + £;sym(p)<52 -I- 0{b^). (2) 

With X — 1, for instance, the Egymip) is what pre- 
dicted by a Hartree-Fock calculation using the Gogny 
effective interaction The parameters Ai{x) and 

Au{x) are Ai{x) = -120.57 + 2Bx/ {a + 1) and Au{x) = 
—95.98 — 2Bx/{<T + 1), respectively. Including also the 
well-known contribution from nucleon kinetic energies, 

i.e., E^^^ip) = (2^/3 _ 1)3^0 (^/^^)2/3 ^ 13.0(p/po)2/3^ 



TABLE L The parameters F (MeV), G , K^ym (MeV), L 
(MeV), and /sTasy (MeV) for different values of x. 



X 


F 


G 


K 


L 


Kasy 


1 


107.232 


1.246 


-270.4 


16.4 


-368.8 





129.981 


1.059 


-88.6 


62.1 


-461.2 


-1 


3.673 


1.569 


94.1 


107.4 


-550.3 


-2 


-38.395 


1.416 


276.3 


153.0 


-641.7 



the density dependence of nuclear symmetry energy is 
shown in Fig. ^for a; = 1, 0, —1 and —2. 

The interaction part of nuclear symmetry energy can 
be parameterized by E'^°^{p) = F{x)p/po + (18.6 — 
F{x)){p/po)^^^\ with F{x) and G{x) given in Table 
12 for X = 1, 0, —1 and —2. Also shown are other 
characteristics of the symmetry energy, including its 
slope L 
9p§ {d^E, 



3po {dEsy^/dp)p^p^ 



and curvature Kgym = 



asy 



at po, as well as the isospin- 
6L I20II of the isobaric 



sym 



dependent part K.^ 
incompressibihty of asymmetric nuclear matter K{5) =^ 
Ko + K^yS^ 0,143. 

For comparisons we have also constructed the follow- 
ing momentum-independent potential (SBKD) that has 
Kq = 200 MeV Jill and exactly the same Esymip) as the 
MDI interaction: 



UsBKDip,S,T) = -356 p/po + 303 (p/po)^/'' 
+ 4t£;p;^(p) + (18.6-F(x)) 
X (Gix) - l)(p/po)«(-)<52. 



(3) 



Isospin diffusion in heavy ion collisions can in principle 
be studied by examining the average isospin asymmetry 
of the projectile-likc residue in the final state. Since re- 
actions at intermediate energies are complicated by pree- 
quilibrium particle emission and production of neutron- 
rich fragments at mid-rapidity, differences of isospin dif- 
fusions in mixed and symmetric systems are usually used 
to minimize these effects jlOl- To study isospin diffusion 
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Sn reactions at = 50 MeV/nucleon and 



an impact parameter of 6 = 6 fm, we thus also consider 
the reaction systems ^^^Sn + ^^^S and ^-^^Sn -I- ^^^Sn and 
^^"^Sn -I- -'^^^Sn at same energy and impact parameter as 
in Ref. . The degree of isospin diffusion in the reaction 
^^''Sn + ^^^Sn is then measured by [53 



Ri 



2X 



i24Sn+ii2Sn 
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(4) 



where X is the average isospin asymmetry (S) of the 
^^''Sn-like residue defined as the composition of nucle- 
ons with local densities higher than po/20 and velocities 
larger than 1/2 the beam velocity in the cm. frame. A 
density cut of po/8 is found to give almost same results. 
In ideal case, the value of Ri ranges between 0.05 and 1 
from complete mixing to full transparency. 
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FIG. 2: (Color online) Time evolutions of Ri and average 
central density for MDI and SBKD interactions with x = —1. 
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FIG. 3: (Color online) Symmetry potential as a function of 
momentum (a) or density (b) with the MDI interaction and 
SBKD interaction (c). 



As an example of present more realistic calculations, 
we show in Fig. |21 the time evolutions of Ri and aver- 
age central density calculated with x = —1 using both 
MDI and SBKD interactions. It is seen that the isospin 
diffusion process occurs mainly from about 30 fm/c to 
80 fm/c corresponding to average central density from 
about 1.2po to O.Spo- However, the value of Ri still 
changes slightly with time until after about 120 fm/c 
when projectile-like and target-like residues are well sep- 
arated. This is partly due to the fact that the isovector 
potential remains appreciable at low density as shown in 
Fig. O where the symmetry potential (C/„ — Up)/2S is 
shown as a function of momentum (panel (a)) or density 
(panel (b)) for the MDI interaction and as a function 
of density for the SBKD interaction (panel (c)). Also, 
evaluating isospin diffusion Ri based on three reaction 
systems, that have different time evolutions for the pro- 
jectile residue as a result of different total energies and 
numbers of nucleons, further contributes to the change of 
Ri at low density. For the two interactions consider here, 
the main difference between the values for Ri appears in 
the expansion phase when densities in the participant re- 
gion are well below po- The experimental data from MSU 
are seen to be reproduced nicely by the MDI interaction 
with X ~ ^1, while the SBKD interaction with x — —1 
leads to a significantly lower value for Ri as the strength 
of the momentum-independent potential is stronger (see 
Fig. |2Jl, which has been shown to enhance the isospin 
diffusion ^lHll^. 

To see how isospin diffusion depends on the density 
dependence of nuclear symmetry energy, we show in Fig. 
0]the final saturated value for 1 — Ri, which measures the 
degree of isospin diffusion, as a function of i^asy for both 
MDI and SBKD interactions. It is obtained by averaging 
the value of 1 — i?^ after 120 fm/c with error bars cor- 
responding to its dispersion, whose magnitude is similar 
to the error band shown in Ref.jS^ for the theoretical 
results from the BUU model. For the SBKD interac- 
tion without momentum dependence, the isospin diffu- 
sion decreases monotonically (i.e., increasing value for 
Ri) with increasing strength of Xasy as the correspond- 



ing isovector potential is mostly positive and decreases 
with increasing stiffness of i^sym (p) in the whole range of 
considered x parameter. The isospin diffusion is reduced 
when the momentum-dependent interaction MDI is used 
as the momentum dependence weakens the strength of 
symmetry potential except for x — —2. As seen in Fig. 
the symmetry potential in the MDI interaction has 
the smallest strength for x = —1 as it is close to zero at 
fc « 1.5 fm~^ and p/po ~ 0.5, and increases again with 
further hardening of the symmetry energy, e.g., x — —2, 
when it becomes largely negative at all momenta and 
densities. The MDI interaction with a; = — 1 thus gives 
the smallest degree of isospin diffusion among the inter- 
actions considered in present study and reproduces the 
MSU data as already shown in Fig. El 
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FIG. 4: (Color online) The degree of isospin diffusion as a 
function of A'asy with the MDI and SBKD interactions. 7 is 
the parameter for fitting the corresponding symmetry energy 
with Esymip) = 31.6(p/po)^- 



The symmetry energy in the MDI interaction with 



~1 is Esymip) = 13.0(p/po)'/^ + 3.7p/po + 
1" w 31.6(p/po)^ °^. It leads to a value of 
-550 MeV for the isospin dependent part of the 
isobaric incompressibility of asymmetric nuclear matter. 
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which should be compared to the pubhshed constraint of 
-566 ± 1350 < ifasy < 139 ± 1617 MeV extracted earher 
from studying giant monopole resonances [l^ . 

It is worthwhile to mention that if the isoscalar part 
of the SBKD potential in Eq. © is replaced with the 
momentum-dependent MDYI potential of Gale et al. 
[33| . which has a similar Jfo as those for the MDI and 
SBKD potentials, the resulting Ri — 0.37 ± 0.07 is much 
closer to that obtained with the MDI = 0.44 ± 0.05) 
than with the SBKD [R, = 0.19 ± 0.06) interaction for 
a; = — 1. The strongly repulsive momentum-dependent 
isoscalar potential thus reduces the effect of isovector po- 
tential on the reaction dynamics. Results from present 
study are therefore not much affected by the uncertainty 
in the momentum dependence of the isovector potential. 

In summary, we have used an isospin- and momentum- 
dependent transport model to study isospin diffusion in 
heavy-ion collisions at intermediate energies. We find 
that the diffusion of isospins happens mostly during ex- 
pansion stage when the density is below normal nuclear 
matter density. The momentum dependence in the nu- 
clear potential plays an important role and affects the 
sensitivity of the degree of isospin diffusion to the density 
dependence of nuclear symmetry energy. Aside from the 
uncertainty due to residual effects of preequilibrium par- 
ticle emission and cluster formation, the present study 
shows that within the context of the mean-field inter- 
actions considered here recent experimental data from 
NSCL/MSU on isospin diffusion is consistent with a nu- 
clear symmetry energy that is nearly linear in density 
at subnormal densities. This leads to a significantly 
constrained value of about —550 MeV for the isospin- 
dependent part of the isobaric incompressibility of isospin 
asymmetric nuclear matter. 
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